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摘要 :CASA 模型 是 研究 宏观 尺度 净 初 级 生产 力 最 常用 的 模型 之 一 ,最 大 光 能 转化 率 是 CASA 模型 的 最 关键 参数 ,但 他 难以 通过 
测量 和 试验 获得 ,只 能 通过 模拟 求 取 。 利 用 CASA 模型 反 演 了 辽河 三 角 洲 河口 湿地 典型 芦苇 群落 的 最 大 光 能 转化 率 , 并 针对 遥 
感 和 气象 数据 的 可 能 误差 对 最 大 光 能 转化 率 的 影响 进行 了 敏感 性 分 析 。 模 拟 结果 表明 :芦苇 群 落 具 有 极 高 的 碳 转 化 能 力 ,最 大 
光 能 转化 率 达 1.667 g CZMJ ,实际 转化 率 达到 0.957 一 1.102 g CZMJ。 敏 感性 分 析 结 果 显 示 最 大 光 能 转化 率 模 拟 值 对 总 辐射 和 
NDVI 表现 出 较 强 的 敏感 性 ,由 总 辐射 误差 带 来 的 最 大 光 能 转化 率 相对 变化 幅度 仅 为 -4.14% 一 4.56% ;模拟 结果 对 NDVI 的 敏 
感性 随 着 NDVI 的 增加 而 降低 ,即便 是 以 误差 30% 考 虑 ,模拟 值 仍然 比较 集中 在 样 点 的 变化 范围 之 内 ,这 些 结果 表明 模拟 的 芦 
苇 最 大 光 能 转化 率 具有 一 定 的 稳定 性 和 可 靠 性 。 
关键 词 :最 大 光 能 转化 率 ;芦苇 ;CASA 模型 ; 净 初 级 生产 力 ; 辽 河 三 角 洲 
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Abstract: Maximum light conversion efficiency ( MLE) is a critical parameter for the Carnegie Ames Stanford Approach 
(CASA) model, which is widely used for modeling net primary productivity ( NPP) globally. However, it is difficult to 
parameterize MLE using experiments and field observation. MLE is fundamental in ecological studies, therefore, modeling 
MLE is of vital importance and significance. The present study determined the MLE of a Phragmites salt marsh in the Liaohe 
River estuarine wetland in China. The main objectives of this study were to: (1) determine the MLE of a Phragmites salt 
marsh; and (2) investigate the sensibility of MLE to environmental factors. Factors included in the CASA model comprised 
the Absorbed Photosynthetically Active Radiation ( APAR) using sunshine duration, which was obtained from the National 
Meteorological Information Center; Fraction of Absorbed Photosynthetically Active Radiation (FAPAR), calculated using 
the Normalized Difference Vegetation Index ( NDVI ) determined from 16-day Moderate Resolution Imaging 


Spectroradiometer ( MODIS) data; temperature and water stress coefficients, calculated using MODIS reflectance data; and 
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meteorological variables including air temperature, precipitation, sunshine duration, air pressure, water vapor pressure, 
wind velocity, and relative humidity. Field observations of the NPP of the Phragmites salt marsh were conducted at 54 
sampling areas with a size of 250 m x 250 m, of which 27 sites were used for modeling, and the remaining for validation. 
The results showed that the Phragmites salt marsh had a relatively high carbon conversion efficiency, with an average MLE 
of 1.667 gC/MJ, varying between 1.112 and 2.611 g C/MJ, which was much higher than the proposed value of 0.389 g C/ 
MJ by Potter. It was even higher than that of broad-leaved, coniferous, and theropencedrymion forests. The simulated MLE 
was sensitive to global solar radiation and NDVI, decreasing with their increase, which was more pronounced at lower 
values. This indicated that it is important to check data quality and increase the data accuracy of global solar radiation. In 
the present study, global solar radiation was estimated using sunshine duration, with an accuracy of 95%. The relative range 
of MLE affected by the error of global solar radiation was from —4.14% to 4.56%. MLE became less sensitive as NDVI 
increased, whereas the simulated values still fell into the MLE range, but the NDVI error increased by 30%. In practical 
applications, the differences in MODIS NDVI data were much smaller, suggesting that our results are universal and could be 
used for other satellite images with different spatial resolutions. The air temperature and precipitation errors had little effect 
on the simulated results, as MLE was not sensitive to them. The results of the sensitivity analysis increased the reliability 
and confidence of the simulated MLE for the Phragmites salt marsh, which is of great significance when studying the carbon 


sink and sequestration potential of Phragmites wetlands in China and other regions globally. 
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River estuarine wetland 


河口 湿地 是 一 类 特别 的 湿地 生态 系统 ,位 于 陆地 海洋 和 河流 生态 系统 的 交界 处 ,是 生产 力 最 高 的 生态 系 
统 之 一 ,可 以 与 农田 生态 系统 相 媲美 ,具有 很 高 的 固 碳 能 力 思 。 单 位 面积 盐 沼 . 红 树 林 湿 地 甚至 比 成 熟 的 热 
带 雨 林 能 封存 高 得 多 的 碳 , 在 全 球 碳 循环 中 扮演 着 重要 角色 。 传 统 的 样 地 观测 是 湿地 生产 力 研 究 的 基本 方 
法 ,站 点 实测 数据 精度 高 ,是 模型 校正 和 验证 不 可 缺少 的 基础 资料 ;但 由 于 湿地 特殊 自然 环境 ,加 大 了 野外 测 
量 的 难度 ; 另 一 方面 , 样 地 观测 在 区 域 净 初 级 生产 力 研 究 中 存在 明显 局 限 性 。 在 区 域 或 全 球 尺 度 上 ,模型 估 犁 
成 为 净 初 级 生产 力 一 项 重要 手段 习 。 其 中 ,以 基于 光 能 转化 率 的 Carnegie Ames Stanford Approach ( CASA ) 模 
型 应 用 较为 广泛 。CASA 模型 能 够 与 RS 和 GIS 技术 相 结合 ,为 区 域 生 产 力 的 研究 提供 高 效 、 精 确 的 研究 
方法 。 

最 大 光 能 转化 率 是 CASA 模型 最 重要 的 参数 之 一 , 它 的 大 小 直接 影响 到 净 初 级 生产 力 的 总 量 。 其 含义 为 
在 理想 情况 下 ,植被 通过 光合 作用 吸收 半 位 光合 有 效 辐射 所 固定 的 最 大 干 物质 总 量 。Raymond 等 认为 在 没有 
气候 和 其 它 因素 限制 情况 下 ,森林 的 最 大 光 能 转化 率 可 达 3.5 g CAMJ' ,而 在 Harvard 森林 研究 样 地 ,最 大 光 
能 转化 率 约 为 0.55 g CAMJ' 。 农 作物 的 最 大 光 能 转化 率 变化 范围 较 窗 , 介 于 1.1 一 1.4 g C/MJ 之 间 ,而 自然 
生态 系统 变化 较 大 六 ,因此 ,Heimann 等 采用 平均 值 1.25 g CZMJ 估算 了 全 球 植被 的 年 兆 初 级 生产 力 !%1 。 不 同 
植被 类 型 之 间 的 最 大 光 能 转化 率 存在 较 大 差异 ,目前 还 无 法 直接 测量 ,只 能 通过 模拟 来 求 取 。Potter 提出 
CASA 模型 时 ,利用 Raich 等 和 McGuire 等 :观测 的 17 组 净 初 级 生产 力 数据 模拟 出 最 大 光 能 利率 为 0.389 
g C/AMJ' 1。 虽然 这 个 值 被 广泛 用 来 估算 区 域 植被 净 初 级 生产 力 '" ,但 最 大 光 能 利率 取 值 偏 低 导致 兆 生 产 力 
被 低估 的 现象 也 被 大 量 报道 ”1 ,这 些 结果 表明 对 最 大 光 能 转化 率 的 本 地 校正 或 者 模拟 对 准确 估算 净 初 级 
生产 力 具 有 十 分 重要 的 作用 。 由 于 受 实测 资料 的 限制 ,相关 研究 仅 有 少量 的 报道 。Running 等 根据 生态 生理 
过 程 模型 BIOME-BGC ,模拟 出 不 同 植被 的 最 大 光 能 转化 率 为 0.389 一 1.259 g CAMJ' SI。 朱文 泉 利用 中 国 林业 
部 的 林业 普查 的 实测 数据 和 空间 分 状 率 为 8 km 探 路 者 数据 ,采用 CASA 模型 模拟 出 中 国 主要 和 森林 植被 的 最 
大 光 能 转化 率 为 0.542 一 0.985 g CZMJI 9 。 由 于 陆地 生态 系统 的 复杂 性 ,对 于 不 同 植被 类 型 ,最 大 光 能 转化 率 
取 值 应 该 不 同 ,森林 生态 系统 的 模拟 值 不 能 简单 应 用 到 湿地 生态 系统 ;另外 ,使 用 的 数据 及 方法 的 差异 也 可 能 
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对 结果 有 一 定 的 影响 。 即 便 不 考虑 这 些 因素 ,迄今 为 止 , 仍 未 发 现 有 关 河 口 湿 地 最 大 光 能 转化 率 模拟 研究 的 
报道 。 

辽河 三 角 洲 河口 湿地 是 我 国 重要 的 河口 湿地 之 一 ,是 亚洲 最 大 的 暖 温 带 滨海 湿地 ,位 于 辽宁 省 西南 部 辽 
河 平原 南端 ,渤海 辽东 湾 的 顶部 (图 1) 。 受 海洋 和 陆地 交互 淡 咸 水 交汇 作用 ,复杂 动力 机 制造 就 了 复杂 多 样 
的 湿地 类 型 和 生态 环境 ,发 育 了 以 芦苇 , 翅 碱 莲 ,茸毛 香 蒲 等 作为 建 群 种 的 湿地 植物 群落 ,以 及 水 稻田 . 虾 蟹 
养殖 池 等 人 工 湿 地 类 型 ,其 中 芦苇 面积 达到 756 km ,是 亚洲 第 一 大 芦苇 湿地 ,已 建成 世界 第 二 大 苇 场 ,修建 
了 完善 的 水 利 排灌 体系 ,改变 了 自然 状态 下 河口 湿地 的 水 文 规律 。 石 油 开发 的 各 种 地 面 工程 造成 了 湿地 景观 
破碎 化 ,严重 破坏 了 湿地 的 原 有 生境 。 在 这 些 人 为 干扰 下 ,辽河 三 角 洲 河口 湿地 的 植被 及 净 生 产 力 格 局 发 生 
明显 的 变化 “i。 但 由 于 种 种 原因 ,我 国 对 湿地 的 研究 监测 开展 较 晚 ,湿地 生态 系统 缺少 长 期 .系统 的 观测 数 
据 。 尽 管 目前 已 经 在 辽河 三 角 洲 开展 了 一 些 净 初级 生产 力 的 研究 工作 "1 ,但 绝 大 部 分 仍 局 限于 有 限 样 点 的 
野外 调查 工作 。RS 和 GIS 技术 为 湿地 净 生 产 力 的 研究 提供 了 高 效 手段 ,但 目前 在 本 区 域 的 应 用 也 还 停留 在 
景观 层面 的 研究 上 "… 。 因 此 ,辽河 三 角 洲 河口 湿地 净 初 级 生产 力 的 遥感 监测 被 列 为 中 国 地 质 调 查 项 目 及 国 
土 资 源 公 益 性 行业 科研 专项 的 重要 内 容 , 而 准确 确定 芦 尘 群落 的 最 大 光 能 转化 率 是 净 初 级 生产 力 研究 的 前 提 
工作 ,对 全 面 了 解 河 口 湿地 碳 汇 功能 、 研 究 湿 地 碳 循环 具有 十 分 重要 的 作用 ,而 且 对 准确 评估 我 国 芦 苇 湿 地 的 
固 碳 功 能 发挥 和 提高 湿地 生态 系统 的 固 碳 潜力 具有 重要 的 意义 。 
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1 辽河 三 角 洲 河口 湿地 的 位 置 及 样 点 的 分 布 
Fig.l Location of the sampling sites in Liaohe river estuarine wetland 
* 背景 为 TM432 波段 合成 


1 数据 与 方法 


1.1 方法 
1.1.1 CASA 模型 

利用 光 能 转化 率 原理 , Monteith 提出 利用 吸收 光合 有 效 辐 射 和 光 能 转化 率 来 估算 陆地 净 初 级 生产 力 的 概 
念 ,并 考虑 了 水 分 .温度 和 养分 胁迫 。 在 此 基础 上 ,Potter 和 Field 建立 了 CASA 模型 ,实现 了 基于 光 能 转化 率 
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原理 的 净 初 级 生产 力 估算 ,其 理论 构架 可 表示 为 下 式 ”1 。 
NPP(%,t) = APAR(x,t) x e(x,t) (1) 
E(x,t)= Togl(x,t) Xx Te2(x,t) x Woe(x,t) XeE (2) 
式 中 ,NPP 为 净 初 级 生产 力 ,APAR 为 吸收 光合 有 效 辐射 ,FPAR 为 光合 有 效 辐射 分 量 ,t 表示 时 间 ,x 表示 空间 
位 置 , se, 为 最 大 光 能 转化 率 ,Te, 和 Te, 为 温度 胁迫 系数 ,We 为 水 分 胁迫 系数 。 
1.1.2 CASA 模型 参数 化 
植被 吸收 的 光合 有 效 辐射 取决 于 太阳 总 辐射 和 植被 对 光合 有 效 辐 射 的 吸收 分 量 , 用 下 列 公 式 表 示 : 
APAR(x,i) = Rs(x,1) xXaXxFPAR(C>，D (3) 
式 中 ,Rs(x, 1) 是 像 元 x 处 的 太阳 总 辐射 , a 表示 植被 所 能 利用 的 太阳 有 效 辐射 占 太阳 总 辐射 的 比例 ,主要 是 
波长 范围 在 0.38 一 0.76 pm 的 可 见 光 , a 一般 约 占 太阳 总 辐射 能 量 的 50% 左 右 。 由 于 太阳 辐射 站 点 数量 少 ， 
常常 需要 利用 采用 Angstrim-Prescott (A-P) 模 型 和 日 照 时 数 来 计算 5 。 
光合 有 效 辐 射 分 量 是 植被 对 光合 有 效 辐射 的 吸收 比例 ,取决 于 植被 类 型 和 植被 覆盖 状况 。 研 究 表明 光合 
有 效 辐射 吸收 分 量 与 归 一 化 植被 指数 (NDVI) 存在 线性 关系 中 。 这 一 关系 可 以 根据 某 一 植被 类 型 NDVI 的 
最 大 值 和 最 小 值 以 及 所 对 应 的 光合 有 效 辐射 吸收 分 量 最 大 值 和 最 小 值 来 确定 : 
(NDVI(%,1) -=-NDVI ，) 
FPAR(x,t) = 2 x (FPAR,. - FPAR,,) +FAPR (4) 
NDVT; »,. — NDVT an 
式 中 ,NDVI, ,和 NDVI ,分别 对 应 第 ;种 植被 类 型 的 NDVI 的 最 大 值 和 最 小 值 。 进 一 步 的 研究 表明 ,光合 
有 效 辐 射 吸收 分 量 与 比值 植被 指数 (RVI) 也 存在 较 好 的 线性 关系 '” ,可 由 公式 表示 : 
Cs (FPAR -FPAR,) + FAPR, (5) 
RVT — RVL,;, i En i 
式 中 ,FPAR,,, 和 FPAR,,, 的 取 值 与 植被 类 型 无 关 , 分 别 为 0. 001 和 0. 95。 研 究 表 明 由 NDVI 所 估算 的 光合 有 
效 辐 射 吸收 分 量 比 实测 值 高 ,而 由 RVI 所 估算 的 光合 有 效 辐射 吸收 分 量 则 低 于 实测 值 ,但 其 误差 小 于 由 
NDVI 所 估算 的 结果 ,考虑 到 这 种 情况 ,Los 将 这 两 种 方法 结合 起 来 , 取 其 平均 值 作为 最 终结 果 '” ,此 时 , 佑 
算 值 与 实测 值 之 间 的 误差 达到 最 小 。 
光 能 转化 率 是 CASA 模型 中 最 关键 的 环节 。Potter 等 认为 在 理想 条 件 下 植被 具有 最 大 光 能 转化 率 ,而 在 
现实 条 件 下 的 光 能 转化 率 受 到 环境 因子 如 气温 水 分 状况 以 及 大 气 水 汽 压 差 等 因素 的 影响 ,因此 利用 温度 、 大 
气 及 水 分 的 可 利用 程度 来 调节 。 其 公式 可 以 表达 为 1. 


FPAR(x,1) 


E(x,t)= Tg (x,t) Xx Tgs(x,t) x We(%,i) x evo (6) 
式 中 ,Teg,(x, t) 和 Tg,(%, i) 为 温度 胁迫 系数 ,分 别 用 以 下 公式 计算 . 
Toei(x) = 0.8 + 0.02Topt(x) - 0.0005Topt (x)? (7) 
1.1814 
Tg,(x) 二 (1+ 02 Topt(#) -10-T(x.0) )(1+ E03 -Topt(#) -10+7(4,0) ) ) (8) 


式 中 ,7 为 一 年 内 NDVI 达到 最 高 时 月 份 的 平均 气温 ,认为 此 温度 为 植被 生长 的 最 适 温度 。 当 某 一 月 平均 温 
度 小 于 或 等 于 -10C 时 ,7 取 0, 认 为 此 时 光合 生产 为 零 。 若 某 月 气温 比 最 适宜 温度 高 10% 或 低 13%C 时 ,该 月 
的 Tg, 值 等 于 月 平均 气温 为 最 适宜 温度 了, 时 的 Tg, 值 的 一 半 。 

We (x, 7) 为 水 分 胁迫 系数 ,反映 了 水 分 条 件 对 光 能 转化 率 的 影响 ,其 取 值 范围 为 0.5 至 1, 由 下 式 计算 . 

We(x,t) =0.5 + 0.5EET(x,1)/PET(x,:) (9) 

式 中 ,EET 为 实际 蔡 散 ,根据 FAO 推荐 的 作物 潜在 燕 散 的 Penman-Monteith 方程 计算 '”。PET 为 潜在 蔡 散 
量 , 根 据 Boucher 提出 的 互补 关系 计算 :2 。 

根据 公式 (1) 和 (2) , 则 最 大 光 能 转化 率 表示 为 : 

因 NPP (x,t) 

™ APAR(x,t) x Te(x,t) x To,(x,t) x We(x,t) 


[34 


(10) 
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公式 中 可 以 看 出 , 当 NPP APAR Te, Te 和 We 为 已 知 时 , 便 可 以 计算 某 一 样 点 上 植被 的 最 大 光 能 转化 率 。 
由 于 模型 估算 往往 不 可 避免 的 要 对 某 些 过 程 进行 简化 ; 另 一 方面 ,不 同 的 数据 误差 也 有 可 能 被 带 入 最 后 的 结 
果 中 ,不 同样 点 所 模拟 的 最 大 光 能 转化 率 可 能 不 同 。 因 此 ,本 文 将 利用 54 个 样 点 数据 中 的 27 个 (图 1) ,根据 
误差 最 小 原则 , 拟 合 出 芦苇 的 最 大 光 能 转化 率 。 
1.1.3 净 初 级 生产 力 分 配 

部 分 净 初 级 生产 力 在 生长 期 分 配 到 根部 ,经 过 累积 形成 地 下 生物 量 。 而 利用 CASA 模型 所 估算 的 净 初 级 
生产 力 包括 地 上 和 地 下 部 分 。 由 于 湿地 的 特殊 困境 ,加 上 芦苇 发 达 的 深 根系 ,对 地 下 净 初 级 生产 力 的 测定 比 
较 困 难 。 本 文 野外 样 地 调查 仅仅 测定 了 六 苇 地 上 部 分 。 为 了 使 CASA 模型 模拟 的 结果 与 实测 生物 量具 有 可 
比 性 ,采用 根 、 茎 、 叶 动态 配置 模式 (Pierre's model) ,根据 可 利用 资源 光 (L) .水 分 ( 丈 ) 和 假 营养 (V) 对 不 同 器 
官 进 行 矶 的 动态 配置 .1 ,将 地 上 部 分 的 净 初 级 生产 力 从 CASA 模型 中 所 模拟 的 结果 分 离 出 来 。 


NPP,, = 3Nr(L/(L + 2min( W,N))) (11) 
NPP = 3Ns(min(W,N)/(2L + min( W,N)) (12) 
NPP = 1 ~ (NPP,,, + NPP,.,) (13) 


式 中 ,Nr 和 Ns 代表 没有 资源 限制 时 , 净 初 级 生产 力 分 配 到 根 和 茎 中 的 比率 ,默认 取 值 0.3。 资 源 可 利用 性 因 
子 工 WN 取 值 介 于 0.1( 资 源 限制 ) 与 1( 资 源 充分 并 容易 获取 ) 之 间 。L 利用 叶 面 积 指数 来 估算 , L =e*” ， 
KK 为 消光 系数 ,默认 取 值 为 0.5。 水 分 因子 下 取决 于 土壤 湿度 和 土壤 质地 。 假 营养 因子 N 利用 湿度 和 气温 来 
计算 ,N= FpxFi, 其 中 fp=PPT/APET,Ft=2[ (7-30)/10] ,PPT 为 降水 量 ,7 为 平均 气温 ,PET 潜在 蒸 散 。 
1.2 数据 收集 与 处 理 

基于 CASA 模型 估算 植被 净 初 级 生产 力 涉及 大 量 的 气象 数据 遥感 数据 以 及 其 它 基础 资料 ,数据 质量 的 
好 坏 直接 关系 到 结果 的 可 靠 性 ,因此 ,在 数据 的 收集 与 处 理 过 程 中 ,要 尽 可 能 的 降低 数据 本 身 所 包含 的 误差 ， 
并 使 数据 具有 统一 的 标准 。 
1.2.1 气象 数据 

气象 数据 包括 2011 年 每 日 的 总 辐射 日 照 时 数 .平均 气温 .最 高 气温 .最低 气 温 大气压、 相对 湿度 。 所 有 
气象 数据 都 从 中 国 气象 科学 数据 共享 网 (http:// cde. cma. gov. cn/) 获取 。 根 据 CASA 模型 估算 区 域 的 净 初 
级 生产 力 , 需 要 基于 网 格 点 的 气象 数据 ,这 要 求 将 气象 站 点 的 气象 要 素 通 过 空间 插值 到 区 域 ,但 研究 区 域内 仪 
有 营口 站 的 气象 数据 ,为 了 减少 插值 的 边缘 误差 ,下 载 了 整个 辽宁 省 内 的 27 个 气象 站 的 数据 ,然后 通过 
ANUCLIM 插值 方法 内 插 出 整个 辽宁 省 网 格 点 气象 数据 ,最 后 利用 研究 区 域 的 边界 剪 切 出 区 域内 的 数据 :1 。 
1.2.2 ”遥感 数据 

遥感 数据 是 进行 区 域 交 初级 生产 力 研 究 不 可 缺少 的 数据 源 。 由 于 MODIS 多 波段 数据 可 以 同时 提供 有 关 
陆地 植被 大 气 .温度 等 特征 信息 ,其 时 间 分 辩 率 高 , 重 覆 周期 短 ,是 开展 区 域 净 初 级 生产 力 研究 非常 重要 的 
资源 。 因 此 ,本文 也 使 用 MODIS 数据 产品 来 模拟 芦苇 群落 的 净 初 级 生产 力 ,具体 包括 2011 年 全 年 16d 合成 
的 MODIS09 反射 率 产品 .2011 年 全 年 16d 合成 MOD13 植被 指数 产品 以 及 2011 年 全 年 16d 合成 MODIS15 叶 
面积 指数 产品 。 所 有 MODIS 数据 均 从 NASA(http://reverb.echo. nasa.gov) 下载 获 取 。 这 些 数据 产品 都 经 过 
了 一 系列 严格 的 质量 控制 和 检验 ,已 经 在 国内 外 广泛 使 用 ,实践 证 明 这 类 产品 的 精确 度 能 够 满足 此 类 研究 工 
作 的 需要 。 
1.2.3 ” 净 初 级 生产 力 测量 数据 

净 初 级 生产 力 测量 采用 收割 法 ,选取 连续 分 布 的 且 长 势 均 质 的 典型 芦苇 群落 ,设置 样 地 54 个 (图 1) ,其 
中 27 个 用 于 模拟 最 大 光 能 转化 率 , 另 一 半 用 于 验证 模拟 结果 。 大 样 地 面积 设 为 230 mx250 m, 以 匹配 MODIS 
像 元 大 小 ,在 每 一 样 地 的 四 角 和 中 心 设置 5 个 2 mx2 m 小 样 方 ,将 样 方 芦 革 齐 地 收割 , 按 器 官 部 位 分 别称 总 
鲜 重 ,并 各 取样 品 一 份 (300 一 500 g 左右 ) ,用 于 测定 不 同 器 官 的 干 物 质 比 率 。 并 收集 样 方 中 所 有 调 落 物 并 测 
定 鲜 重 ,同样 也 取样 品 一 份 ,测算 调 落 物 的 干 物质 比率 。 将 样品 带 回 实验 室 在 80%C 下 烘 干 48 bh 至 恒 重 并 称 
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重 , 然 后 根据 不 同 器 官 的 干 物质 比率 及 其 总 鲜 重 ,推算 其 样 方 内 总 干 重 , 以 5 个 小 样 方 的 平均 值 计 算 样 地 单位 
面积 生物 量 ,最 后 乘 以 碳 转换 系 数 (0.45) ,换算 成 单位 面积 上 植被 地 上 部 分 碳 含量 。 各 样 点 的 调查 具体 情况 
见 表 1。 


表 1 模拟 样 点 的 最 大 光 能 转化 率 


Table 1 The simulated maximum light conversion efficiency of Phragmites salt marsh 


模拟 样 点 编号 


最 大 光 能 转化 率 /(g CZMJ ) 


Saplings 经 度 /( .) 纬度 人 植被 柳 章 度 观测 NPPZ(kg C/m?) Masimum light 
sites No Longitude Latitude Vegetation coverage Observated NPP Bonverston eficieney 
1 121.64584 40.91238 95% 0.604 1.757 
2 121.68842 40.94095 90% 0.493 1.291 
3 121.68153 40.96292 85% 0.771 1.930 
4 121.73628 40.95840 90% 0.424 1.122 
5 121.71923 41.00459 80% 0.490 1.607 
6 121.75887 41.00171 80% 0.664 1.712 
7 121.80473 41.02351 60% 0.402 1.158 
8 121.72478 41.02634 80% 1.090 2.297 
9 121.75261 41.04469 75% 0.435 1.177 
10 121.82411 41.04255 70% 0.483 1.193 
11 121.78935 41.07868 70% 0.463 1.112 
12 121.75301 41.11547 80% 0.588 1.667 
13 121.74062 41.16428 80% 0.633 1.738 
14 121.65253 41.23127 80% 0.540 1.397 
15 121.70103 41.26432 70% 0.406 1.182 
16 121.73387 41.20943 70% 0.534 1.386 
17 121.80019 41.20012 70% 0.872 2.013 
18 121.80817 41.15675 80% 0.779 2.037 
19 121.88129 40.91645 80% 0.978 2.610 
20 121.86287 40.93732 75% 0.879 2.256 
21 121.93085 40.94338 70% 0.707 1.804 
22 121.93343 40.96785 70% 0.868 2.353 
23 121.89491 40.96758 75% 0.679 1.780 
24 121.89799 40.98792 60% 0.590 1.572 
25 121.84760 41.20659 60% 0.482 1.208 
26 121.84058 41.25302 75% 0.496 1.330 
27 121.80091 41.24509 80% 0.589 2.031 


2 结果 与 讨论 


2.1 模拟 结 

表 1 列 出 了 模拟 的 芦苇 最 大 光 能 转化 率 , 样 点 模拟 值 为 1.112 一 2.611 g CZMJ (1.656+0.43) ( 表 1) ,模拟 
值 平均 值 为 1.667 g CZMJ。Potter 提出 CASA 模型 时 ,植被 最 大 光 能 转化 率 为 0.389 g CZMJ” ,明显 低 于 本 
区 域 芦苇 群落 的 最 大 光 能 转化 率 。 在 广东 易 湖 山 和 黑石 顶 的 观测 也 表明 常 绿 益 叶 林 的 实际 光 能 转化 率 分 别 
要 大 于 0.605 g CZMJ 和 0.515 g CAMJ' ”1 ,因此 , 彭 少 麟 采用 最 大 光 利 用 率 为 1.25 g CZMJ 估算 了 广东 省 植被 
的 净 初 级 生产 力 '”1 ,但 依然 要 低 于 芦苇 群 落 的 最 大 光 能 转化 率 。 采 用 与 本 文 相 同 的 方法 ,朱文 泉 利 用 净 初 级 
生产 力 的 实测 数据 ,根据 误差 最 小 的 原则 模拟 出 中 国 主要 森林 植被 类 型 的 最 大 光 能 转化 率 为 0.542 一 0.985 g 
CZMJI20 。 而 Running 根据 生态 生理 过 程 模 型 BIOME-BGC 模拟 的 不 同 的 植被 类 型 的 最 大 光 能 转化 率 为 
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0.389 一 1.259 g C/AMJ' 1 ,其 结果 也 都 明显 低 于 本 区 域 芦 苇 群 落 。 
利用 模拟 的 最 大 光 能 转化 率 , 我 们 计算 了 芦苇 群落 
的 实际 光 能 转化 率 , 样 点 模拟 平均 值 为 0.957 一 1.102 g 
CZMJ ,要 高 于 我 国 植被 的 年 平均 光 能 转化 率 0.196 一 
0.701 g CZMJI5 ,甚至 要 高 于 常 绿 益 叶 林 、 常 绿 针 叶 林 、 
热带 常 绿 阔 叶 林 等 森林 植被 "1 。 图 2 为 芦苇 群落 光 能 
转化 率 的 时 间 动 态 , 表 现 出 明显 的 单 峰 曲线 季节 变化 模 
式 ; 从 年 初 到 4 月 上 旬 ,由 于 低温 天 气 的 限制 , 光 能 转化 
率 趋 近 于 0; 当 温度 升 高 到 适宜 植被 生命 活动 时 ,植被 TO 
开始 生长 发 育 , 光 能 转化 率 开 始 上 升 ; 随 着 生长 季 的 持 I Day of year ld 
续 , 温 度 和 降水 的 增加 有 利于 叶片 的 光合 作用 , 光 能 转 a 
化 率 不 断 提高 , 叶 面 积 不 断 增加 ,， 黎 盖 度 逐渐 增 大 ,在 7 。 Fig2 Temporal vartion of light coversion efficiency of 
月 下 旬 本 地 区 的 温度 和 降水 分 别 达到 最 大 值 时 , 光 能 转 phragmites salt marsh 
> 化 率 也 达到 最 高 值 1.524 g CZMJ; 在 植被 长 势 达 到 最 茂 
2 盛 后 , 随 着 温度 和 降水 的 减少 ,加 上 大 部 分 绿叶 开始 变 黄 或 调 落 ,叶绿素 含量 降低 , 光 能 转化 率 也 开始 降低 ;在 
一 ”生长 季 结 束 后 ,低温 寒冷 天 气 限制 了 植被 的 生长 , 光 能 转化 率 趋 近 于 0。 
OO 河口 湿地 位 于 陆地 海洋 和 河流 生态 系统 的 交界 处 ,是 生产 力 最 高 的 生态 系统 之 一 局。 乔治 亚 海岸 互 花 
米 草 的 净 初 级 生产 力 达 到 7620 g C/m*'* ,我 国 长 江口 典型 芦 革 湿 地 净 初 级 生产 力也 可 达 1100 一 2400 g C7/ 
nm。 这 除了 与 河口 湿地 所 处 的 特殊 环境 有 关外 ,还 可 能 取决 于 湿地 植被 高 的 光 能 转化 率 。 本 文 的 结果 表 
1 明 芦 革 群 落 具 有 很 高 的 最 大 和 实际 光 能 转化 率 , 甚 至 要 高 于 Running “|, 朱文 泉 ““ 模 拟 的 阔 叶 林 、 针 叶 林 以 
© 及 针 阔 混交 林 等 森林 植被 ,这 可 能 是 由 于 植被 的 生长 都 需要 充足 的 水 分 和 营养 条 件 ,而 芦苇 沼泽 湿地 恰恰 处 
于 陆地 与 水 体 的 过 渡 地 带 , 它 同时 兼 具 丰 富 的 陆 生 和 水 生动 植物 资源 ,在 长 期 的 自然 选择 和 进化 过 程 中 形成 
了 其 它 任何 生态 系统 都 无 法 比拟 的 基因 库 ; 别 一 方面 ,沼泽 湿地 的 泥 碳 土壤 经 过 长 期 累积 ,为 沼泽 植被 的 生长 
提供 了 丰富 的 营养 元 素 ;加 上 河口 湿地 海洋 与 河流 的 水 文 交 互 作用 有 利于 植被 光合 作用 对 碳 的 转化 ,已 有 的 
研究 也 表明 河口 湿地 高 的 生产 力 与 水 文 的 季节 性 波动 有 着 密切 的 相关 性 ' 引 。 这 些 特殊 的 水 文 和 土壤 条 件 为 
河口 湿地 芦苇 群落 对 碳 的 转化 提供 了 有 利 的 条 件 , 形 成 了 芦苇 群落 高 的 光 能 转化 率 。 
() ”3.2 模拟 结果 验证 1.6 
对 模拟 结果 进行 验证 是 进一步 分 析 和 应 用 的 基础 ， 
由 于 在 目前 还 不 能 通过 试验 和 设备 来 测量 植被 的 最 大 
光 能 转化 率 ,无 法 通过 实测 值 来 检验 模拟 的 效果 。 最 大 
光 能 转化 率真 是 CASA 模型 的 重要 参数 ,因此 ,利用 54 
个 野外 测量 数据 中 的 另外 27 个 ,并 利用 模拟 的 芦苇 最 
大 光 能 转化 率 来 计算 对 应 样 点 的 地 上 净 初 级 生产 力 , 通 
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过 验证 地 上 净 初 级 生产 力 来 间接 验证 模拟 的 最 大 光 能 0 0 

关 a 、 | 人 直人 直 扑 0 0.4 0.8 1.2 
转化 率 的 可 靠 性 。 图 3 为 CASA 模型 模拟 的 芦 革 地 上 et 
净 初 级 生产 力 与 实测 值 散 点 图 。 总 体 上 看 ,CASA 模型 


模拟 结果 与 观测 值 具体 较 好 的 一 致 性 。 模 拟 值 与 实测 ei 

值 之 间 的 相关 系数 为 0.63( P<0.001) | 达到 极 显著 性 水 J Scatter plots of aboveground NPP estimation vs. 
observation 

平 。 最 大 相对 均 万 根 误差 (RRMSE) 为 40% ,平均 相对 ， 原始 结果 是 指 采 用 最 大 光 能 转化 率 0.389 g CZMJ 所 计算 的 

均 方 根 误差 为 24%。27 个 验证 点 的 平均 模拟 值 和 观测 。 NPP ,模拟 结果 是 指 采用 1.667g CZMJ 所 计算 的 NPP 

值 较为 接近 ,分 别 为 619.9 g C/m* 和 627.3 g C/m*, 两 
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者 仅 相 差 7.5 g CMm: 。 这 间接 的 表明 本 文 所 模拟 的 芦 著 最 大 光 能 转化 率 具 有 较 高 的 可 靠 性 。 而 当 最 大 光 能 
转化 率 采用 0.389 g CZMJ 时 ,模拟 的 地 上 净 初 级 生产 力 为 112.87 一 184.59 g CZMJ (146.95+13.96) (图 3) , 明 
显 低估 了 本 区 域 芦苇 群落 的 固 碳 能 力 。 中 国 芦 苇 资 源 丰 富 ,分 布 集中 而 广泛 ,全 国 范围 内 有 14 个 主要 芦苇 产 
区 ,芦苇 湿地 面积 达 130 万 hm200 。 按 照 本 文 的 最 大 光 能 转化 率 计算 , 则 芦苇 的 年 固 碳 能 力 可 达到 10.04 一 
15.14 Tg C(13.55+1.25) ;而 最 大 光 能 转化 率 采用 0.389 g C/AMJ 时 ,模拟 的 芦苇 年 净 初 级 生产 力 仅 为 2.34 一 
3.53 Tg C(3.16+0.29) , 则 芦苇 湿地 的 固 碳 能 力 将 会 被 低估 75% 。 
3.3 模拟 的 最 大 光 能 转化 率 对 主要 因素 的 敏感 性 

利用 CASA 模型 来 模拟 植被 的 最 大 光 能 转化 率 , 需要 先 利 用 遥感 数据 和 气象 数据 计算 APAR ,FPAR ,ze 
和 Tg, 因此 ,这 些 数 据 误差 也 会 被 带 和 人 最 后 的 模拟 结果 中 。APAR 的 数据 误差 主要 来 源 于 利用 日 照 时 数 对 太 
阳 总 辐射 的 估算 ,FPAR 的 误差 则 主要 来 源 于 NDVI, We 和 Teg 则 主要 取决 于 气温 和 降水 。 将 这 些 主要 的 影响 
因素 以 2011 年 的 实际 值 为 基准 ,在 +40% 的 变化 范围 内 ,以 5% 的 间距 变化 ,采用 单 因 素 分 析 误 差 对 模拟 结 
的 影响 。 由 于 植被 的 NDVI 的 值 介 于 0 一 1 之 间 , 当 变化 超过 30% 时 ,出 现 NDVI 大 于 1 的 情况 ,因此 NDVI 的 
变化 范围 设 为 :30%。 

图 4 结果 为 不 同 影 响 因 素 值 域 下 模拟 的 芦苇 最 大 光 能 转化 率 和 最 大 光 能 转化 率 对 其 影响 因素 的 一 阶 导 
数 。 从 图 中 可 以 看 出 ,不 同 的 因素 对 模拟 结果 的 影响 不 同 。 当 总 辐射 和 NDVI 在 +40% 的 范围 内 变化 时 ,模拟 
的 芦苇 最 大 光 能 转化 率 分 别 为 1.239 一 2.568 g CZMJ(1.752+0.409) 和 1.313 一 2.276 g CAMJ(1.716+0.308) , 模 
拟 值 随 着 总 辐射 和 NDVI 的 增加 而 呈现 下 降 的 趋势 (图 4) ,并 表现 出 较 强 的 敏感 性 ,其 一 阶 导数 分 别 为 
-1.351 一 -0.642(-0.917+0.23) 和 -1.218 一 -0.708(-0.924+0.17) 。 在 低 值 区 域 ,这 种 敏感 性 表现 的 更 加 明 


显 , 并 虽 着 总 辐射 和 NDVI 的 增加 而 呈现 下 降 的 趋势 (图 4) 。 尽 管 如 此 ,模拟 值 仍然 落 在 样 点 的 变化 范围 之 
内 (1.112 一 2.611 g CZMJ) 。 本 文 的 结果 表明 在 高 纬度 及 总 辐射 较 低 的 地 区 ,控制 总 辐射 的 精度 显得 格外 重 
要 ,因为 利用 CASA 模型 模拟 的 植被 最 大 光 能 转化 率 可 能 具有 更 大 不 确定 性 。 但 由 于 观测 设备 和 观测 方法 的 
限制 ,我 国 太 阳 辐 射 观测 站 稀少 , 常 需要 利用 日 照 时 数 来 估算 。 在 本 研究 区 域 也 没有 太阳 总 辐射 观测 站 ,只 能 
利用 日 照 时 数 来 计算 '”。 即 便 是 整个 辽宁 省 也 仅 有 朝阳 、 大 连 和 沈阳 具有 辐射 观测 数据 ,因此 ,利用 这 3 个 
点 的 数据 ,采用 交叉 验证 的 方法 来 检验 所 估算 的 太阳 总 辐射 的 精度 。 检 验 结 果 显 示 估 算 的 太阳 总 辐射 具有 很 
高 的 精度 (图 5) ,3 个 样 点 的 相对 均 方 根 误差 均 小 于 7% ,总 体 精度 达到 95% , 按 +5% 的 误差 估计 , 则 模拟 的 芦 
化 最 大 光 能 转化 率 介 于 1.598 一 1.743 g CZMJ 之 间 ,绝对 变化 幅度 为 -0.069 一 0.076 g CZMJ ,相对 变化 幅度 仅 
为 -4.14% 一 4.56%。 这 表明 本 文 模拟 的 结果 具有 一 定 的 稳定 性 和 可 靠 性 。 

模拟 的 芦苇 最 大 光 能 转化 率 对 NDVI 也 具有 和 较 强 的 敏感 性 ,NDVI+30% 的 变化 引起 的 模拟 结果 绝对 变化 
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图 4 模拟 的 最 大 光 能 转化 率 对 主要 因素 的 敏感 性 


Fig.4 Sensibility of the simulated maximum light conversion efficiency to affect factors 
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图 5 太阳 总 辐射 的 观测 值 与 预测 值 
Fig.S Scatter plots of global solar radiation estimation vs. observation 
RRMSE : 相对 均 方 根 误差 Relaitve root mean square error 
幅度 为 -0.35 一 0.609 g CZMJ, 相对 变化 幅度 为 -21. 1.00 ee 
CA rd m 一 一 m 
2% 一 36.5% ,这 也 证 明 空 间 分 辨 率 和 植被 覆盖 度 对 模拟 
的 结果 有 重要 影响 。 随 着 植被 覆盖 度 的 增加 ,这 种 敏感 ba 
性 有 降低 的 趋势 ,如 图 4 所 示 , 一 阶 导数 随 着 NDVI 的 对 vs 
增加 而 呈 上 升 ( 负 值 ) 的 趋势 。NDYVI 的 不 确定 性 是 个 有 
十 分 复杂 的 问题 ,有 来 自 不 同 的 遥感 平台 、 传 感 锅 、 遥 感 0.25 
方式 的 差异 的 影响 ,即便 是 不 考虑 这 些 因 素 , 云 .大气 成 
、 总 世 0.00 
分 ,地面 背景 观测 角度 等 因素 也 都 对 NDVI 有 着 重要 1 49 97 145 193 241 289 337 
、 2 a 、 寺 序 Day of year/d 
的 影响 。 刘 良 明 等 指出 MODIS NDVI 的 不 确定 性 约 为 
12% , 随 着 NDVI 值 的 增加 ,这 种 不 确定 性 还 会 减少 :52 ; 图 6 两 种 分 辨 率 的 芦苇 MODIS NDVI 


而 Trishchenko 等 的 研究 表明 MODIS NDVI 数据 的 误差 Fig.6 Temporal vartion of MODIS NDVI with the resolution of 
高 达 30%'”。 以 最 大 误差 30% 考虑 , 则 模拟 结果 为 人 ma 
之 内 (1.112 一 2.611 g C/AMJ)。 别 一 方面 ,由 于 MODIS 
NDVI 被 广泛 作为 区 域 植 被 及 净 初 级 生产 力 的 研究 的 首选 数据 源 。 相 同 的 数据 源 使 得 这 些 因素 引起 的 NDVI 
的 不 确定 性 具有 相似 性 ,使 最 终结 果 之 间 差 异 将 会 大 大 减 小 ,结果 具有 可 比 性 。 即 便 是 使 用 不 同 尺度 或 者 分 
辨 率 的 MODIS NDVI 数据 (图 6) ,其 结果 也 十 分 相近 。 从 图 6 中 可 以 看 出 ,250 m 和 1000 m 分 辩 率 的 MODIS 
NDVI 几乎 重合 ,差异 约 为 1% ,模拟 值 变化 仅仅 为 0.016 g CAMJ。 这 是 因为 MODIS NDVI 的 合成 都 是 利用 分 
辩 率 为 250 m 的 第 一 波段 和 第 二 波段 反射 率 为 基础 ,优先 选择 近 星 下 点 无 云 像 元 , 尽 可 能 减 小 云 暗影 大气 
气 溶胶 的 影响 .应 用 BRDF 模式 将 角度 订正 到 一 致 的 标准 化 的 视角 和 太阳 大 顶 角 的 条 件 下 ,保证 合成 资料 的 
质量 和 一 致 性 659 。 因 此 ,本 文 基于 250 m 分 辨 率 的 MODIS NDVI 数据 的 结果 仍然 可 以 应 用 到 1000 m 分 辩 率 
尺度 下 。 

模拟 的 最 大 光 能 转化 率 随 着 气温 的 增加 呈 上 升 的 趋势 ,但 随 着 降水 量 的 增加 呈 下 降 的 趋势 ,两 者 的 一 阶 
导数 都 较 小 ,并 且 保 持 比 较 恒定 的 值 ( 图 4) ,表明 模拟 结果 对 温度 和 降水 量 的 敏感 性 较 小 而 且 比 较 稳 定 。 模 
拟 值 分 别 为 1.547 一 1.801 g C/AMJ (1.669+0.08) 和 1.586 一 1.774 g C/AMJ (1.672+0.06) 之 间 , 绝 对 变化 幅度 分 
别 为 -0.12 一 0.134 g CAMJ 和 -0.08 一 0.107 g C/AMJ。 而 实际 上 ,温度 和 降水 的 观测 简单 而 且 精 度 高 ,其 误差 远 
远 要 小 于 50%。 因 此 ,温度 和 降水 误差 对 模 结 果 的 影响 更 加 微弱 。 
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准确 确定 最 大 光 能 转化 率 是 利用 光 能 利用 率 模型 研究 区 域 植被 净 初 级 生产 力 的 基础 。 本 文 利 用 遥感 

术 结 合 地 面 气象 观测 和 生产 力 野 外 监测 ,利用 CASA 模型 反 演 了 辽河 三 角 洲 河口 湿地 典型 芦苇 群落 的 最 大 光 
能 转化 率 。 芦 苇 群 落 具 有 很 高 的 碳 转 化 能 力 , 最 大 光 能 转化 率 达 到 1.667 g CZMJ ,实际 光 能 转化 率 达 到 
0.957 一 1.102 g C/AMJ。 模 拟 最 大 光 能 转化 率 对 总 辐射 和 NDVI 表现 出 较 强 的 敏感 性 ,在 低 值 区 域 ,这 种 敏感 
性 表现 的 更 加 明显 ;因此 ,在 总 辐射 较 低 的 地 区 ,控制 总 辐射 的 精度 显得 格外 重要 。 本 文 总 辐射 模拟 精度 达到 
95% , 按 +5% 的 误差 估计 , 则 模拟 值 相对 变化 幅度 仅 为 -4.14% 一 4.56% ,具有 一 定 的 稳定 性 和 可 靠 性 。 模 拟 结 
果 对 NDVI 的 敏感 性 随 着 NDVI 的 增加 而 降低 ,即便 是 以 误差 30% 考 虑 , 则 模拟 值 仍 然 比 较 集 中 在 样 点 的 变 
化 范围 之 内 。 而 在 实际 应 用 中 ,由 于 MODIS NDVI 被 广泛 用 于 区 域 净 初 级 生产 力 的 研究 ,不 同 分 辩 率 的 
MODIS NDVI 差异 较 小 ,模拟 结果 也 更 加 稳定 。 模 拟 芦 革 最 大 光 能 转化 率 对 温度 和 降水 量 的 敏感 性 较 小 而 且 
比较 稳定 ,因此 ,温度 和 降水 误差 对 模拟 结果 影响 较 小 。 
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